Analysis of pion elliptic flows and HBT interferometry in a granular quark-gluon 

plasma droplet model 
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In many simulations of high-energy heavy-ion colUsions on an event-by-event analysis, it is known 
that the initial energy density distribution in the transverse plane is highly fluctuating. Subsequent 
longitudinal expansion will lead to many longitudinal tubes of quark-gluon plasma which have 
tendencies to break up into many spherical droplets because of sausage instabilities. We are therefore 
motivated to use a model of quark-gluon plasma granular droplets that evolve hydrodynamically 
to investigate pion elliptic flows and Hanbury-Brown-Twiss interferometry. We flnd that the data 
of pion transverse momentum spectra, elliptic flows, and HBT radii in ^sjvjv ~ 200 GeV An 4- 
Au collisions at RHIC can be described well by an expanding source of granular droplets with an 
anisotropic velocity distribution. 

PACS numbers: 25.75.-q, 25.75.Nq, 25.75.Gz 
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I. INTRODUCTION 

Recently, there has been much progress in our un- 
derstanding of the process of nucleus-nucleus collisions 
at RHIC 11 0, H S S H- Hydrodynamical calcu- 
lations agree well with the RHIC V2 data of the el- 
liptic flow at low transverse momentum < 2 GeV 
0, H, En Ol- However, thejy cannot predict the 
saturation of V2 at higher p^, H, 0, liol llll and the 
RHIC HBT puzzle of i?out/i?side « 1 Ullll El El- 
The HBT puzzle is contrary to many earlier theoreti- 
cal expectations [S El El El US iSl • Various mod- 
els have been pu t forth to explain the HBT puzzle 
111 1110 IS IH 111 IS 111 113 111 111 111 IP . There 
have also been many attempts to provide a consistent 
explanation for both the elliptic flow and HBT measure- 
ments at RHIC §^^^^M,M^- 

In Ref. [23, a granular particle-emitting source of 
quark-gluon plasma (QGP) droplets evolving hydrody- 
namically was put forth to explain the RHIC HBT puz- 
zle. The suggestion was based on the observation that in 
the hydrodynamical model [TgL |4^ , the particle emission 
time scales with the initial radius of the droplet. Par- 
ticles will be emitted earlier if the radius of the droplet 
is smaller, as in a source of many droplets. An earlier 
emission time will lead to a smaller extracted HBT radius 
i?out) while the extracted HBT radius i?sido is determined 
by the scale of the distribution of the droplets. As a re- 
sult, the value of i?out can be close to i?side for a granular 
quark-gluon plasma source |28|- Further suggestions of 
using the single-event intensity interferometry to search 
for the signature for granular structures have also been 
presented lisLli^ . 

Motivated by the successes of our previous analysis of 
the HBT puzzle using a granular droplet model, we exam- 
ine in this paper the theoretical basis for the possible oc- 
currence of granular structure in the evolution of a quark- 



gluon plasma. In addition, we wish to refine the granular 
source model of Ref. by considering more reasonable 
anisotropic velocity distributions of the droplets instead 
of a constant radial velocity assumed in [2^1 . We would 
like to investigate the pion elliptic flows as well as HBT 
radii as functions of the pion transverse momentum. The 
agreement both of the V2 and HBT radii with experi- 
mental data concurrently will gives strong constrains for 
the granular model and its parameters and will provide 
a useful insight into the initial state and the evolution 
of the particle-emitting source produced in high-energy 
nucleus-nucleus collisions. 



II. GRANULAR INSTABILITY IN THE 
EVOLUTION OF THE QUARK-GLUON PLASMA 

Based on the recent results of high-energy heavy-ion 
collisions at RHIC, the matter (presumably QGP) pro- 
duced in the collisions may be a strongly-coupled medium 
with a very high energy density H H ^ lilllM; They 
are thermalized within about 1 fm/c ^|,11 H H H 0- 
It is of interest to study its subsequent evolution and see 
whether granular structures may play a role in the space- 
time development of the presumed quark-gluon plasma 
matter. 

Although a granular structure was suggested earlier as 
the signature of a first-order phase transition [i^, the 
occurrence of granular structure may not be limited to 
the occurrence of a first-order phase transition. There 
arc additional effects which may lead to the dynamical 
formation of granular droplets. 

Because the QGP is a strongly interacting dense 
medium, a surface tension arises at a boundary due to 
the presence of a strong interaction in the dense phase on 
one side of the boundary and the absence (or weakening) 
of the strong interaction with no density (or diminishing 



density) on the other side |M|47llMlM|50l|5ll|5a|53|. 

This imbalance of the forces acting on different parts at 
the boundary leads to the surface tension, while the de- 
tail profile of the boundary may depend on the nature 
and the order of phase transition. Due to the presence 
of this surface tension, there may be instability against 
surface shape changes and bulk density oscillations. 



A. Transverse Density Fluctuations and Surface 
Granular Instability 

The quark gluon plasma that may be produced in the 
laboratory is spatially bounded by various boundaries. 
Furthermore, in many simulations of the heavy-ion colli- 
sions on an evcnt-by-evcnt basis, the initial transverse 
energy density is far from being uniform. It exhibits 
large transverse density fluctuations with a large peak- 
to- valley ratio 22j Jd, .Mj ■ A large number of transverse 
density domains are clearly visible in the initial trans- 
verse density distribution of the produced matter in Fig. 
21 of or Fig. 4 of 0. Because of these highly 
fluctuating initial transverse density distribution form- 
ing "lumps" in different regions, the probability for the 
occurrence of granular structure may be enhanced. 

We can present an approximate scenario to discuss the 
evolution of matter with a large transverse density fluc- 
tuations, in a central collision of two heavy and equal 
nuclei of radius R. Because of Lorentz contraction, the 
longitudinal length i?/7 of the produced matter is much 
smaller than the transverse length, which is of the order 
R. Initially as depicted in Fig. 1(a), density fluctuations 
in the transverse plane manifest as a number of trans- 
verse lumps characterized by a radius rdn , which one can 
also see from the results of and |4l|. From such an 
initial data, we can consider the evolution according to 
Landau hydrodynamics 'ssj , which gives a good descrip- 
tion of the longitudinal distribution of the produced par- 
ticles [5^ ■ In such a description, because the initial 
longitudinal length is much smaller than the transverse 
length, the hydrodynamical force per unit volume in the 
longitudinal direction, —dP/dz, is much greater than the 
hydrodynamical force per unit volume in the transverse 
direction, ^dP/dp. The subsequent expansion in the 
longitudinal direction proceeds much faster than in the 
transverse direction, as pointed out by Landau and Be- 
lenkii 55]. The expanding matter develops into tubes 
of radius r^o [Fig. 1(6)] and the density of matter in the 
tube decreases as a function of the proper time. 

As is well known, a long tube of matter is unstable 
against "sausage" perturbations that tend to make the 
tube narrow in some regions and thicker in other regions. 
The sausage instability leads to the break-up of the long 
tube of matter into many approximately equal spherical 
droplets, as in the break-up of toroidal liquid droplets 

We can examine the sausage instability for a cylinder 
with a length L and a radius r^o, with L » rdo- We 



(a) Initial configuration 



(b) Shortly afterwards 



(c) Later afterwards 
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FIG. 1: Schematic description of the matter distribution un- 
der a rapid longitudinal expansion in the ^-direction, for a 
highly fluctuating initial transverse density distribution [Fig. 
la]. The system evolves into many longitudinal tubes at Fig. 
lb and later at Fig. Ic. Sausage instabilities lead to the pro- 
duction of droplets along the longitudinal tubes as depicted 
in Fig. Ic. 



consider "sausage" perturbations a„ about the cylinder 
of matter, 



rdiz) = fdo<l+ ^ ' 

k n— odd 



an sin(n7rz/L) 



(1) 



In order to have a positive slope drd{z)/dz at z = and 
a negative slope aX z — L, we need to limit ourselves to 
perturbations with odd values of n and a„ > 0. By the 
condition of volume conservation, the quantity f^o is 



4 ^ 

n— odd 



^^ n— odd 

The surface energy of the cyhnder is then 



n— odd 



Es = cj27rrdo^|l - 



n— odd 




The cylinder is unstable against perturbations of order n 
if (PEs/da^ is negative. We can evaluate d'^Es/daf^ and 
obtain 



d^Es/dai = cr27rrrfoL{-l/2+ (2/n7r)2}. 



(4) 



Thus, d^Es/da^-^ is negative for all odd integer values 
of n. As an unrestrained growth of the perturbation 
of order n leads the formation of {n + l)/2 droplets, 
we find that the cylinder is unstable against perturba- 
tions leading to the formation of 1, 2, 3, ... approximately 
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equal droplets. Thus, when the length L of the tubes are 
stretched, sausage (granular) instability will develop and 
the long tube will fragment into approximately equal- 
size droplets, as depicted in Fig. 2c. Large density fluc- 
tuations in the transverse direction, together with the 
surface tension effects, favor the formation of granular 
droplets. 

In our simplified description, we have ignored the ef- 
fects of the perturbation on the variation of the trans- 
verse surfaces at the two ends of the cylinder. A refine- 
ment on the description of the end surfaces will not mod- 
ify significantly the question of the instability of the cylin- 
drical matter, for cylinders with L substantially greater 
than Tdo- 

In conclusion, we find from the above analysis that the 
occurrence of granular droplets may be more common 
than previously thought. It is of great interest to explore 
the subsequent dynamics in a granular model for the de- 
scription of various observables measured in high-energy 
heavy-ion collisions. This is particularly relevant because 
a granular droplet model can explain the occurrence of 
early emission of particles and a small value of _Rout ■ 



B. Bulk density oscillations due to surface tension 
interactions 

To carry out the analysis of the effects of surface ten- 
sion on the stability of density oscillations, we consider 
the following simple model of the surface tension, as rep- 
resented by a finite-range interaction. The equation of 
hydrodynamics without a finite-range interaction is js^ 



duj 



dp 

dx'- 



UiU 



dp 



(5) 



where e and p are the energy density and pressure. In 
this hydrodynamical equation, interactions between con- 
stituents leading to bulk properties of the plasma have 
already been included into the characterization of the 
pressure p and the energy density e. What remains is the 
effective residual finite-range interaction ^(r — r') that 
can lead to the surface tension. In the presence of this 
finite-range residual interaction, the above can be gener- 
alized to be [63 



[e(r) + p{r)]u^{v) 



dui{r) dp{r) 



- / dr 



w w 



Viv{r 



dx^ 
-r') 



dp{r) 
dx'^ 

(6) 



where W is the average energy per particle which may 
be temperature-dependent, and e{r)/W gives the number 
density n{r) of the constituents at r. 

We choose to represent the effect of the surface ten- 
sion in terms of an effective residual interaction of the 
following form as used in li^. 



/ir e 



fj.r 



(7) 



where > and 1/fi are the strength and the range of 
the residual interaction that leads to the surface tension. 
As a residual interaction, v{r) has been chosen such that 
the contribution of the interaction to the total energy of 
a system with a uniform density is zero and is positive 
and proportional to the surface area in the surface region. 

To carry out the stability analysis, we can envisage 
that by solving the hydrodynamical equation ^ without 
the surface interaction, one obtains a comparatively slow- 
varying density distribution eo{r,t) with a local average 
Co = (eo(r, i)). Then following the derivations of the 
equation of sound waves [H^, the perturbation Se{r,t) = 
e(r, t) — eo(r, t) obeys the wave equation 



d^de{r,t) 



cl\/^6e{r,t)+Vr ■ J dr'[{Se{r, t)e{r',t) 
+ e(r,^)(^e(r',^))]Vr^;(r-r')/W^^ (8) 



where Cg is the speed of sound. We consider density per- 
turbations of the type 



^e(r, t) = e 



ik-r— ZLjf 



(9) 



Representing the spatial variation of eo(r,f) by its aver- 
age eo, we obtain from Eq. ||HJ) 

= c2/(k)k2 , (10) 

for points in the interior of the medium. Here, /(k) is 



/(k) = 1 



(k2 + /i2)2 ' 



(11) 



where the second term of /(k) arises from the surface 
tension interaction v(r), and a is the granular stability 
number defined as 



I^Vcf 



(12) 



Equations and Hll|l indicates that when f{k) be- 
comes negative, the system is unstable against a density 
perturbation of wave number k. We plot in Fig. 2 the 
function /(fc) as a function of k/n. As one observes, the 
minimum of /(fc) is located at k/fi = 1 for which the 
characteristic wave length of the oscillation 1/fc is equal 
to the range of the interaction The function f{k) 

becomes negative when a > 1 and it is therefore appro- 
priate to call a the granular stability number. When the 
granular stability number a is greater than unity, the sys- 
tem is unstable against perturbations with wave numbers 
k in the interval fc_ < k < k^, with k± given by 



|2a - 1 ± V(2a- 1)^ - l| 



(13) 



We note from Eq. (|12|l that the condition a > 1 for 
the occurrence of granular instability happens when the 
surface tension interaction strength is large. It will 
be useful to relate to the surface tension coefficient 
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relationship 



0.5 1 1.5 2 2.5 3 3.5 4 

k/|i 



FIG. 2: /(k) curves for different values of a. The expansion 
is unstable in the shadowy region. 



a. For such a purpose, we study the energy of a system 
with a numbers density profile given by 



n(r) = nout + 



1 + e^/"^ ' 



(14) 



where no(> riout) is the number density in the interior 
dense region of a surface and riout is the number density 
in the outer dilute region. The energy of the system due 
to the residual interaction of Eq. Q is 

E=^J dpdzn{r) J dp' dz' n{r')v{r - r'). (15) 

We can introduce V{z) as 

V{z) ^ j dp'dz'n{r')v{v -r'). (16) 

After integration the p' coordinates, we find 

/oo 
dz'n{z') 
-oo 

X |e-^l"-^'l -//|z- z'|e-^l^-^'l| . (17) 

From Eqs. (|15|I - H17|) . the surface tension coefficient is 
therefore given by 

a = dE/dp 

/oo /.oo 
dzn{z) I dz'n{z') 
oo J-oo 

X |e-^l^-^'l - fi\z - z'|e-''l^-^'l } , (18) 

which gives a relation between a and ar for a general 
density profile. For the simple case of a sharp bound- 
ary with n{z) — nQ9{—z), the above equation gives the 



riQirar 
2^3 



(19) 



The surface tension coefficient increases with ar/ and 
the square of the number density. As no varies with 
temperature as while ar and fi are relatively slowly 
varying functions of T, the surface tension varies with 
temperature approximately as T^. When the surface is 
described by a sharp surface profile, the granular stabil- 
ity number a is related to the surface tension coefficient 
a by 



2/i(T 



(20) 



In quenched QCD lattice gauge calculations, the sur- 
face tension coefficient (j{Tc) has been calculated for 
T = Tc, and of T^J/T;? was found to range from 0.016 
to 0.092 US El m m. If we use such a value of the 
surface tension coefficient to evaluate the granular stabil- 
ity number a using Eq. H12|l . the quantity a will have the 
value from 0.024 to 0.140 and the system is stable against 
density oscillations of the type discussed here. However, 
the medium considered by the quenched QCD approxi- 
mation consists of gluons without dynamical quarks. As 
is well known, the presence of dynamical quarks has im- 
portant influences on many thermodynamical properties 
of the plasma. The phase transition temperature is al- 
tered from 269 MeV for quenched QCD without dynam- 
ical quarks to 154 MeV for full QCD with 3 flavors of 
quarks 61], whereas the surface tension coefficient a de- 
pends on the temperature approximately as T^. The 
number density of a quark-gluon plasma with three fla- 
vors of quarks is about two times the number density 
of the plasma without dynamical quarks, while the sur- 
face tension coefficient varies as the square of the number 
density. The surface tension coefficient a calculated in 
quenched QCD at T = Tc may be quite different from 
that calculated in full QCD at higher temperatures. It 
will be of interest to evaluate the surface tension coeffi- 
cient in full QCD and study its temperature dependence 
in order to find out whether the quark-gluon plasma is 
stable against bulk density oscillations with wave lengths 
close the the range of the interaction discussed here. 

The above discussions was carried out in the linearized 
perturbative theory of slow hydrodynamical motion. We 
envisage that the initial QGP matter is highly com- 
pressed in the longitudinal direction and the subsequent 
longitudinal expansion will be rather rapid and non- 
linear in the density changes. In this respect, it is of 
interest to note that in previous calculations in full hy- 
drodynamics for the rapid expansion of a nuclear system 
with a finite range interaction, the one-dimensional frag- 
mentation of the density into one-dimensional "lumps' 
occur, as shown in the left panel of Fig. 2 of Ref. |62. If 
these hydrodynamical calculations can be a useful ana- 
logue, the analogous occurrence of granular droplets in 
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the bulk expansion of the quark-gluon plasma may be 
possible. 



III. MODEL OF GRANULAR QGP DROPLETS 

Based on the above discussions, we would like to use a 
model of granular quark-gluon plasma droplets to exam- 
ine the dynamical evolution of the system. We assume 
that the system consists of Nd droplets of radius ini- 
tially distributed in a short cylinder of length 2TZz along 
the beam direction (z direction) as in Fig. 2c, with an 
initial transverse spatial distribution up to a radius TZt-, 



dP,, 



2np dp dz 



[1 - exp(-pVA7^?)] 

xB{nt- p)e{n,-\z\), (21) 



where p = \J + and z are the coordinates of the 
center of a droplet. The parameter ATZt describes a shell- 
type radial distribution which may arise from the dynam- 
ics of the transverse expansion. A naive blast-wave-type 
expansion with a constant radial velocity would leave a 
void along the central transverse axis at p = 0. Instead 
of a constant transverse velocity, we shall find later that 
the best parameters that describe the experimental data 
suggests a transverse velocity profile {Pd)T ^ p^'^^ . From 
the equation of continuity, the density would change ap- 
proximately by 



dn 
It 



-n--^p{/3d)^ 



Pdp 



oc 



-1.42- 



,0.58 ■ 



(22) 



Thus the density decreases more rapidly near p ^ and 
there is a tendency to form a shell-type radial distribu- 
tion. The shell-type distribution of the droplets may also 
arise from the shadowing of detected particles originat- 
ing from the central region of p ~ 0. The HBT radii 
as a function of the average pion transverse momentum 
of a pion pair turns out to be slightly sensitive to this 
parameter of ATZf 

Because of the early thermalization and the anisotropic 
pressure gradient, the droplets will acquire anisotropic 
initial velocities. We assume that the velocity of a droplet 
depends on the initial coordinates of the droplet center, 
iri,r2,r3) = {x,y,z), in the form 



1,2, and 3, (23) 



i(^d)i = a»sign(r,)( — 



where a.i is describes the magnitude of the anisotropic 
expansion, sign(ri) denotes the sign of r^, and hi (6^ , bz) 
are the exponential power parameter that describes the 
variation of the velocity with r.^ . 

For each single droplet, we use relativistic hydrody- 
namics and the equation of state of the entropy den- 
sity ll^ II^. Is^ . Im| to describe its subsequently evolu- 
tion |28l |4J|. The evolution of the granular source of 
many droplets can then be obtained by superposing all 
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FIG. 3: Negative pion transverse momentum spectrum for 
jyj < 0.5. The open circles give the PHENIX data for 
\/snn = 200 GeV Au -I- Au collisions with minimum bias 
[65| . Various curves give theoretical predictions with differ- 
ent parameters. 



of the evolutions of individual droplets |28l l44j . In our 
calculations, the transition temperature and the transi- 
tion temperature width are taken to be Tc = 165 MeV 
and AT = 0.05Tc d IM El El • The initial energy den- 
sity of the droplets is taken to be cq = 3.75TcSc, which 
is about two times of the density of quark matter at Tc 
0, E3- The pions are emitted out of the surfaces of 
droplets at the freeze-out temperature Tf = O.OSTc, with 
momenta obeying the Bose-Einstein distribution in the 
local frame at the temperature Tf. 

In our model calculations there are five velocity pa- 
rameters, {ax, Qy, Uz, = bx = by, and bz}, and four 
size parameters of the granular source, {TZt, Ti-z, ATZt, 
and rd}. We find that the pion transverse momentum 
spectra and the elliptic flow V2 are sensitive to the ve- 
locity parameters and insensitive to the size parameters. 
While our HBT results are sensitive both to the velocity 
parameters and the size parameters. We first determine 
the velocity parameters by comparing the pion transverse 
momentum spectra and the elliptic flow V2 of our model 
calculations with experimental data. Then, we determine 
the source size parameters by comparing our results of 
HBT radii with experimental data. 

Figure 3 shows the PHENIX tt" transverse momentum 
spectrum in y^sJm — 200 GeV Au + Au collisions with 
minimum bias 65] and the transverse momentum spec- 
trum obtained in our model calculations. The theoretical 
results have been normalized by matching the distribu- 
tion to the first bin. In our calculations, we use a cut 
of particle rapidity \y\ < 0.5 in order to compare with 
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experimental data, which have been normahzed to unit 
rapidity |65| . Instead of the velocity magnitude parame- 
ters of ax and ay, we can equivalently use their average 
= {ax + o,y)/2 and their difference Aa^ = ax — ay. 
The transverse momentum spectra depend on a^ and is 
independent of Aa^ . They also depend on aj, 6^, and bz- 
We shall show that the elliptic flow V2 is very sensitive 
to AflT, and h^. Good agreement with experimen- 
tal data of transverse momentum spectra and the elliptic 
flow can be obtained concurrently by taking b^ = 0.42, 
hz = 0.03 and Aot = 0.10. After b^ and bz are fixed, 
the parameters a^, and can be obtained by compar- 
ing the transverse momentum spectrum of the granular 
sources with experimental data ||65|- We finally deter- 
mine Obj. — 0.365 and az — 0.850 for our granular source. 



IV. RESULTS OF THE ELLIPTIC FLOW 



0.3 
0.25 
0.2 



0.15 - 



0.1 

0.05 




b,=0.03, 

bj=0.42, Aa-r=0.10 
b-,=0.03, Aa^=0.06 
bj=1.00, Aaj=0.13 

b^=0.42. b-,=0.42, 
Aa-,=0.13 





(a.,=0.365, a,=0.850) 
• PHENIX DATA 
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The transverse momentum distribution of particles can 
be represented in the form |,6& . .67. .68i] 



dN 

2TTdp^ 



1 



2 Vn cos(n(/ 



(24) 



is the transverse momentum of 



where Pj, = 

the particle, is its azimuthal angle with respect to 
the reaction plane, and the harmonic coefhcients, ?;„ are 
anisotropy parameters. The elliptic flow is defined as the 
second harmonic coefficient V2, which describes the ec- 
centricity of the particle distribution in the momentum 
space. We choose the direction of x-axis in the reaction 
plane and the direction of y-axis out of the reaction plane. 
We can express V2 as 



V2 



(cos 20) 



P^ 



Pt 



(25) 



Figure 4 shows the PHENIX V2 data as a function of pt 
for charged pions in Au + Au collisions at ^/smn = 200 
GeV 10]. Also shown are the theoretical V2 results in 
our granular droplet model. We use the same cut for 
particle pseudo-rapidity, jryj < 0.35, in theoretical calcu- 
lations, as in the experiment set-up (lOj . In the region 
Pj, < 2 GeV, the V2 coefficient for a granular source in- 
creases with p as in previous hydrodynamical calcula- 
tions m m El III IZI Izi- However, in the high 
region with p^ > 2.5 GeV, we find that the V2 co- 
efficient for a granular source with small bz decreases 
with which is different from previous hydrodynam- 
ical calculations H S E IH HI HI Iz3 El- In the 
droplet-rest frame, the average transverse momentum of 
particles in the x and y directions, (p'x) and (p'y), are the 
same. In the laboratory frame, a elliptic flow arises from 
the anisotropic collective transverse boost of the droplet. 
For a very small bz , the longitudinal velocity of droplets 
is almost constant in the source while the magnitude of 
the longitudinal boost velocity az is very large (0.85) in 



FIG. 4: Pion elliptic flow V2 as a function of particle trans- 
verse momentum for the granular sources and the experimen- 
tal data of charged pions of PHENIX fiot . 



our model. In this case, the longitudinal component of 
droplet velocity is important and it does not boost the 
transverse momenta of the particles. The transverse ve- 
locity of the droplet is much smaller than the longitudinal 
velocity of the droplet. Its boost effect becomes weaker 
for the particles with higher transverse momentum p' 
(higher p^). On the other extreme for larger bz, the 
longitudinal velocity of droplet is small for the droplets 
with small longitudinal coordinates. The transverse com- 
ponents of droplet velocity are dominant in the small z 
region, which lead to a production of the elliptic flow and 
the V2 coefficient increasing with transverse momentum 
|7l| , even at higher p^, as shown in the dashed-dot curve 
in Fig. 4. 

After and have been fixed at the values 0.365 
and 0.850 as determined by the pion transverse momen- 
tum spectra, we examine the pattern of V2 as a func- 
tion of transverse momentum px as we vary 6^ , bz, and 
Afly. In the region of < 2 GeV, the curvature of 
V2 is sensitive to 6^. For bz = 0.03, the curvature for 
b^ = 0.03 is too large (see the dashed line in Fig. 4) 
and the curvature for 6^ — 1.00 is too small (see the 
dot hue in Fig. 4). For bz ^ bj, ^ 0.42, the curve 
of V2 rises continuously and deviates from experimental 
data for > 1.5 GeV. We find that the set of param- 



eters {bz,b^,a^} — {0.03,0.42,0.10} give results consis- 
tent with experimental V2 data, as indicated by the solid 
theoretical curves in Fig. 4. The fact that 6^ is much 
greater than bz while is substantially greater than a-p 
indicates that the dynamical behavior in the transverse 
expansion and longitudinal expansion are very different. 

Recent experimental data of the elliptic flow of pion 
show that V2 remains to have the value of 0.1 at very 
high p^ l?^ . The processes of particle production in the 
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high region are dominated by parton and hard-probe 
processes and is beyond the thermal emission model we 
consider here. 



V. RESULTS OF HBT INTERFEROMETRY 

The two-particle Bose-Einstein correlation function is 
defined as the ratio of the two-particle momentum dis- 
tribution P{pi,P2) relative to the the product of the 
single-particle momentum distribution P{pi)P{p2). For 
a chaotic pion-emitting source, P(pi) {i = 1,2), and 
P{pi,P2) can be expressed as [TSf 
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(26) 
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FIG. 5: Two-pion correlation functions for granular source. 



(27) 



where A{pi,Xi) is the magnitude of the amplitude for 
emitting a pion with 4-momentum pi = (p^,Ei) in the 
laboratory frame at Xi and is given by the Bose-Einstein 
distribution with freeze-out temperature Tf in the lo- 
cal rest frame of the source point, ^[pi^p^] Xi, X2) is 
the two-pion wave function. Neglecting the absorption 
of the emitted pions by other droplets, (^{pi,p2] Xi, X2) 
is simply 

^{pi,P2;Xi,X2) 
1 



^/2 



A(pi,Xi)A(p2,^2)e*''^-^^+''^^-^^ 



+ A{p^,X2)A{p2,X^)e'P'■''-+'P-■''' 



(28) 



Using the components of "out" , "side" , and "long" 
[tgL u% of the relative momentum of the two pions, 
q = |pi — P2I, as variables, we can construct the 
correlation function C((7out, gsidc, giong) from P{pi,p2) 
and P(pi)P{p2) by summing over pi and p2 for each 
((Joutjgsidcgiong) bin |2^]. The HBT radii i?out, -Rsido, 
and i?iong can then be extracted by fitting the calculated 
correlation function C{qout, 9sido, ^long) with the following 
parametrized correlation function 

C((?out, gsidc, giong) = 1 + A e-«o„tfl.Lt-9L.«L.-9f„„,fl.f„„, . 

(29) 

Figure 5 shows the theoretical two-pion correlation 
functions for the granular source. The top figures give the 
results for a average pion transverse momentum of a pion 
pair, Kt, less than 400 MeV/c, and the bottom figures 
for Kt > 400 MeV/c. The velocity parameters of the 
granular source are dj. = 0.365, Aa^ = 0.10, — 0.850, 
= 0.42, and bz = 0.03, which are determined by the 
transverse momentum spectra and the elliptic flow V2 as 
discussed above. The size parameters of the granular 
sources are taken to be TZt = 8.8 fm, ATZt — 3.5 fm, 
TZz = 7.0 fm, and rd — 1.3 fm. The number of droplet 
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FIG. 6: Twojjion HBT radii obtained by the PHENIX Col- 
laboration |l4l ] and the STAR Collaboration compared 
with the theoretical results calculated in the granular droplet 
model with ATZt = and 3.5 fm. 



Nd is taken to be 40 in our calculations. One observes 
noticeable differences in the lower Kt region as compared 
to the higher Kt region. 

The left panels of figure 6 give the fitted two-pion HBT 
radii for the granular source as a function of Kt- The 
symbols of circle and down-triangle are for ATZt = 0.35 
fm and ATZt = 0, respectively. The parameter ATZt of 
Eq. (|21|l describes the radial distribution of the centers 
of the droplets, with ATZt = for a uniform distribution 
and ATZt = 0.35 fm for a shell-type radial distribution. 
The experimental PHENIX results [H, and STAR re- 
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suits Q are shown on the right panels. The curve gives 
the theoretical results for ATZt = 0.35 fm. For our theo- 
retical HBT calculations, we use a cut for particle pseudo- 
rapidity region \r]\ < 0.35, the same as in the PHENIX 
experiments jl^ • We find that if we increase the param- 
eter rd, the HBT radii i?out and i?iong will increase. And 
if we increase the parameter ATZt , the variation of HBT 
radii -Rout and i?sidc with Kt will become steep. The 
HBT results of granular source for ATZt = 3.5 fm agree 
quite well with experimental data, although the case with 
ATZ^ = also give almost as good an agreement. 

VI. SUMMARY AND CONCLUSION 

The expansion of the dense matter produced in high- 
energy heavy-ion collisions may favor the production of 
granular droplets. The strong interaction associated with 
the dense matter will lead naturally to the presence of a 
surface tension on a boundary. 

In an event-by-event basis, the initial transverse den- 
sity distribution of matter in central nucleus- nucleus 
collisions have been known to be highly fluctuating 
[2£, .41j This large spatial fluctuations may facil- 

itate the formation of granular droplets. We envisage 
that as the expansion in the longitudinal direction pro- 
ceeds much more rapidly than in the transverse direction, 
these initial transverse density fluctuations will lead to 
the formation of many longitudinal tubes of matter with 
a length much greater than its radius. Due to the surface 
tension effects, these tubes are unstable against granular 
sausage instability and will tend to break up into approx- 
imately equal-size spherical droplets. 

Motivated by the possibility for the occurrence of gran- 
ular droplets and our previous successful HBT analysis in 
terms of granular QGP droplets, we refine our previous 
granular QGP droplet model to include initial anisotropic 
velocities due to the anisotropic pressure gradients at an 
earlier stage of the collisions. We investigate concurrently 



the data of (1) the pion transverse momentum spectra, 
(2) the elliptic flow V2 as a function of , and (3) various 
HBT radii as a function of Kr^ , for Au -I- Au collisions at 
y/^WN = 200 GeV at RHIC. We find that the collection 
of different pieces of experimental data can be described 
well by a granular source model. Direct confirmation 
of the granular droplet configuration may need to await 
more experimental work on single-event HBT correla- 
tions and the fluctuation between correlation functions 
of the single- and mixed-events, as suggested previously 

Ellil. 

The extracted parameters from our analysis provide 
useful information on the nature of the transverse and 
longitudinal expansion. The elliptic flow V2 is due to 
the anisotropic initial dynamical conditions in different 
transverse directions. Its decrease at high is sensi- 
tive to the exponential power parameter of the droplet 
velocity bz in the longitudinal direction. The value of 
bz is much smaller than that of 6^ , indicating different 
dynamical behaviors in longitudinal and transverse direc- 
tions. Our HBT results further indicate that the granu- 
lar source has a small lifetime and a possible shell-type 
spatial distribution. Further studies on the origin of the 
difference in the transverse and longitudinal dynamics as 
well as the spatial distribution in the transverse direction 
will be of great interest. 
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